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We propose a novel dynamical method for beating decoherence and dissipation in open quantum 
systems. We demonstrate the possibility of filtering out the effects of unwanted (not necessarily 
known) system-environment interactions and show that the noise-suppression procedure can be 
combined with the capability of retaining control over the effective dynamical evolution of the open 
quantum system. Implications for quantum information processing are discussed. 
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All real world quantum systems interact with their sur- 
rounding environment to a greater or lesser extent. Such 
systems are said to be open. No matter how weak the 
coupling that prevents the system from being isolated, 
the evolution of an open quantum system is eventually 
plagued by nonunitary features like decoherence and dis- 
sipation Quantum decoherence, in particular, is a 
purely quantum-mechanical effect whereby the system 
loses its ability to exhibit coherent behavior by getting 
entangled with the ambient degrees of freedom. Deco- 
herence stands as a serious obstacle common to all appli- 
cations relying on the capability of maintaining and ex- 
ploiting quantum coherence. These encompass quantum 
state engineering Q , quantum intcrfcrometry || , macro- 
scopic quantum mechanics [Q and, notably, the whole 
emerging field of quantum information processing || . 

Recently, considerable effort has been devoted to de- 
sign strategies able to counteract the effects of environ- 
mental couplings in open-system evolutions. In partic- 
ular, the theory of quantum error correction has been 
developed to meet this challenge From a physical 
point of view, the general underlying question can be 
stated in terms of attaining quantum noise control. Un- 
like the closed-system limit, where the active manipu- 
lation of unitary dynamics is currently realized to be a 
problem of quantum control theory j?J , the possibility of 
directly exploiting control techniques to influence open- 
system properties has not been fully explored yet. Exist- 
ing approaches mainly rely on feedback (or closed-loop) 
control configurations In fact, conventional quan- 
tum error correction protocols can be regarded, in their 
essence, as a form of quantum feedback control imple- 
mented on a redundant physical system. 

In this Letter, we formulate a model for decoupling 
a generic open quantum system from any environmen- 
tal interaction through simpler so-called open-loop con- 
trol techniques. We show that the resulting descrip- 
tion not only provides a comprehensive framework for 



decoherence-suppression schemes as first identified in M 
and subsequently implemented by various authors un- 
der specific assumptions ppj , but, in contrast to previ- 
ous proposals, it also points out a general criterion for 
engineering effective open-system evolutions that are, in 
principle, immune to noise and decoherence. More pre- 
cisely, we find that one can effectively control the system 
to undergo a wide range of dynamical behavior while still 
eliminating the effects of the environment. The allowable 
dynamics are generated by a subgroup of the possible 
system tranformations. This has potentially important 
consequences for quantum control and quantum compu- 
tation, in that it can be regarded as a strategy for per- 
forming fault-tolerant control. Even though the effects of 
the environment make it impossible to retain control over 
arbitrary unitary evolutions of a quantum system, an ef- 
fective dynamics can be still reliably constructed over a 
restricted set of transformations. 

The starting point of the method consists in recog- 
nizing that no relaxation process can take place instan- 
taneously. Accordingly, one should be able to inter- 
fere with the associated dynamics by inducing motions 
into the systems, which are faster than the shortest 
time scale accessible to the reservoir degrees of freedom. 
The usage of tailored time-dependent perturbations as a 
tool to improve system performances has a long history 
within high-resolution Nuclear Magnetic Resonance spec- 
troscopy, where versatile decoupling techniques are avail- 
able to manipulate the overall spin Hamiltonian . De- 
spite this enlighting similarity, the construction of analo- 
gous procedures for open quantum systems faces an im- 
portant conceptual difference, for we assume that any 
decoupling action can be exerted only on the system vari- 
ables, the environment being contributed by a huge num- 
ber of uncontrollable quantum degrees of freedom. 

We consider a quantum system S coupled to an arbi- 
trary bath B, which together form a closed system de- 
fined on the Hilbert space Ti. = TLs ® 'Hb, T~(-s an d Hb 
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denoting S and B Hilbert spaces respectively. The over- 
all Hamiltonian can be written in a concise form as 

H = H s ®l B + t s ®H B + H SB =^2S a ®B a , (1) 



where 1 is the identity operator and the bath operators 
B a are supposed to be linearly independent. Being H 
Hermitian, the linear space spanned by the system op- 
erators S a is a self-adjoint subspace in the vector space 
B(Hs) of bounded operators acting on Hs- We shall 
assume that the unwanted noise-inducing Hamiltonian 
Hsb is only contributed by a finite subset of open-system 
couplings. We call interaction space Tg C B(TCs) the 
corresponding finite dimensional subspace. The second 
ingredient we introduce is the control algebra, Cg, which 
is generated by the repertoire of Hamiltonians we can 
turn on for S to implement decoupling. We allow for 
possibilities where Is ^ Cs- 

If Ptot(O) = ps(0)<8>ps(0) is the initial state over H, the 
open-system evolution of 5* is the coarse-grained dynam- 
ics p s {0) h-> Ps(t)= Tr B {ptot(t)}, Tr B denoting partial 
trace over Hb 0- The relaxation dynamics for ps(t), 
which involves a combination of quantum decoherence 
and dissipation mechanisms depending on the nature of 
the coupling operators, may display a complicate time de- 
pendence. In the simplest case, the off-diagonal matrix 
elements of ps(t) behave like exp(— t/r re i), r re i indicat- 
ing the time scale for significant departure from unitarity 
and irreversible loss of quantum coherence. 

In order to protect the evolution of S against the effect 
of the interaction Hsb, we start by seeking a perturba- 
tion Hi(t) E Cs to be added to Ho as a suitable decou- 
pling interaction, H(t) = Hq + H 1 (t) ® 1 b • We restrict 
here to situations where the control field is cyclic, i.e., as- 
sociated to a decoupling operator U\(t) that is periodic 
over some cycle time T c > 0: 



Ux{t) = Texp 



1 J duHi(u) 



Ut(t + T c 



(2) 



In the interaction representation associated with H\(t), 
defined by ptot (t) = Ui(t)ptot(t)U\(t) on H, time evolu- 
tion is ruled by a transformed time- varying Hamiltonian, 



H{t) = u\{t)H U x {t) = [ul^SMt) 



B a . (3) 



Since U\(T C ) = I5, the evolution in the original 
Schrodinger representation can be easily derived. One 
can prove that after TV cycles, T/v = NT C , 



N 



-iHT N 



(4) 



where the motion of the system under the time- 
dependent field Hi(t) has been replaced by a strobo- 
scope development under an effective so-called average 
Hamiltonian H The calculation of H is performed 
on the basis of a standard Magnus expansion of the 



time-ordered exponential defining the cycle propagator 
U to t(T c )=exp(-iHT c ), 

Texp/-* f TC duH(u)\ =e -il» W +H w +...]T c ^ (5) 



where the various contributions collect terms of equal or- 
der in the transformed Hamiltonian. In particular, 



ff (0) = -L / C duH(u). 

1 c Jo 



(6) 



We shall say that fcth-order decoupling is achieved if 
the control field Hi(t) can be devised so that contribu- 
tions mixing S and B degrees of freedom are no longer 

present in H^ and the first nonvanishing correction 

arises from H , k > 1 . Owing to the fact that the cycle 
time T c enters the Magnus series as a controllable ex- 
pansion parameter, we examine the limit of fast control, 
where the lowest-order terms are expected to provide an 
accurate description (first-order decoupling). Formally, 
for a finite evolution time T, this requires considering 
T c = T/N in the limit as N — > 00. The Magnus series 
defining evolution over a single cycle converges for suffi- 
ciently large N whenever H (r) = 0(T c r ) = 0{1/N 2 ) for 
r > 2. As a result, in the limit of arbitrarily fast control, 
contributions higher than zeroth-order are negligible in 
(H) and we can focus on the problem of designing the 

effective Hamiltonian H . 

We now show that the time average defining H can 
be made identical to a group-theoretical averaging pro- 
cedure [H. Since we can apply any Hamiltonian in the 
control algebra, full control of the decoupling propaga- 
tor Ui(t) is available over the associated set of unitary 
transformations. Let Q be any finite group of unitary 
operators that generates Cs, Q = {gj}, j = 0, . . . , \Q\ — 1, 
\Q\ = ord(Q). Then the map 

S»S = Uc{S) = ^- g\S gj , SeB(H s ), (7) 

is the unique operation projecting on the so-called cen- 
tralizer of Q. Equivalently, since averaged operators S 
commute with every gj, they belong to the so-called 
commutant C of the control algebra. Notice that C is 
closed under commutation. The map (0) is implemented 
through a simple piecewise constant decoupling operator: 



Ui{t)= gj , jAt<t<(j + l)At, 



(8) 



corresponding to a partition of the cycle time T c into \Q\ 
intervals of equal length At = T C /\Q\. Then, by (Q), 



i7 (0) =n c (-ffo) = ^5 Q ®6 ct , 



(9) 



which, by virtue of the quantum operation (Q), displays 
well-defined symmetry properties. 
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The decoupling prescription (JSj) requires the capabil- 
ity of instantaneously changing the evolution operator 
from gj to gj + i over successive subintervals, implying ar- 
bitrarily strong "kicks" of control field. Such impulsive 
full-power control configurations correspond to so-called 
quantum bang-bang controls as introduced in ||. In fact, 
this method can be seen to provide an explicit control 
implementation of an abstract unitary symmetrization 
procedure recently proposed by Zanardi |l3| . 

Two different situations arise depending on the knowl- 
edge available on the system-bath interaction Hsb- Let 
us first suppose that no knowledge is assumed, in which 
case the environmental coupling is completely arbitrary 
and Is = B{TLs). We can prove the following: 

Theorem. Let S be a finite dimensional system and 
let the interaction with the environment be arbitrary, 
Is = B(TLs). Then, in the limit of arbitrarily fast control 
rate, the evolution of observables in the control algebra 
can be suppressed arbitrarily well: 

lim Tr s {A Ps {T = NT C )} = Tr s {Ap s (0)} . (10) 

TV— >oo 

If, in addition, the control algebra is maximal, Cs — 
B{TLs), then complete first-order decoupling is achiev- 
able through system manipulations alone: 

lim p s {T = NT C ) = p s (0) . (11) 

N— >oo 

Proof. Let A 6 C S - The first statement follows by Eq. 
(||) with Hamiltonian (^|) and the fact that L4,iS a ] = 
for every a. If Cs consist of all operators, then the corn- 
mutant C only contains c-numbers, S a — X a Is in (|^) 
and the result follows. □ 

The group able to average system operators into the 
commutant of Cs can be chosen as a set of linearly in- 
dependent unitary operators realizing a so-called unitary 
error basis on TLs |Q- Such subgroups always exist for 
finite dimension. Within the class of piecewise constant 
decoupling sequences as considered above, it can then be 
shown that at least dim(Cg) = \Q\ steps are needed in 
a cycle to attain decoupling. Since the effective system 
evolution is completely quenched by the decoupling pro- 
cedure, we call this configuration maximal averaging. 

When some knowledge is available on the coupling 
Hsb, this information can be exploited to engineer 
shorter decoupling sequences fulfilling specific symmetry 
constraints. For a given (known) interaction space Is, 
the goal is to devise a control algebra able to selectively 
averaging out Is, lie (Is) — 0) while leaving invariant 
the sector in the operator space of the system contain- 
ing some useful dynamics. This requires that the error- 
inducing and the desired system operators in ([j]) trans- 
form according to different irreducible representations of 
Q. By using (|9|), the effective open-system evolution over 
time T is then governed by 

lim p s (T = NT C ) = e- m * T p 8 (0) e+ m * T , (12) 

JV— >oo 



the projected Hamiltonian Hs — Hc(Hs) only consist- 
ing of operators belonging to C. Equivalently, the set of 
operators in the commutant ofCs determines the reliable 
manipulations left available for effective system control. 
For a given Is, the identification of a minimal group Q 
able to produce decoupling is nontrivial. We provide an 
illustrative example. 

Let us discuss a -ftT-qubits dissipative quantum register 
The maximum possible complexity of error gener- 
ation arises in the presence of so-called total decoher- 
ence, whereby combined errors can occur to any num- 
ber of qubits. In this case, I s = B(H S ) ^ (C 2 )® K 
and maximal averaging is demanded to decouple the 
register from quantum noise. Since an error basis on 
q2®k j g g enera ted via the tensor product of the stan- 
dard Pauli bit / sign- flip error basis [jl4j , a possible choice 
is Q = {1 s ,(t^}^ k , a = x,y,z,i = l,...,K. Thus, a 
number of \Q\ = A K steps is required for minimal length 
decoupling sequences. A more efficient averaging is pos- 
sible if the relevant register-bath coupling is known to be 

linear in single-qubit operators, Hsb = J2i Q ® 
Under that condition, which is met for both independent 
decoherence (dim(Is) = 3-JQ and collective decoherence 
(dim(Zs) = 3), selective averaging suffices to decouple 
from errors. One can show that the tensor power of the 

Pauli group, Q = {t s ,®f =1 o- { a ] }, a = x,y,z, \Q\ = 4, 
represents a minimal choice. In terms of the control field 
Hi(t), decoupling is then enacted by cycling the qubits in 
the register through sequences of collective 7r-pulses along 
two axes, e.g., At — n x — At — iT- z — At — 'K- x — At — TT- Zl 
At = T c /4. In the special case of a purely decohering cou- 
pling, Hsb — J2i °z ® B* \ the scheme can be further 
simplified by taking Q = {Is, ®f=iO~x \G\ = 2, which 
reproduces the situation analyzed in [|9| for K = 1. In 
the decoupling regime, the register effectively behaves as 
a noiseless quantum memory, which is an essential ingre- 
dient for various quantum cryptographic and communi- 
cation schemes [§,^5). In addition, one can now perform 
in a fault-tolerant way any logical operation belonging to 
the commutant of the tensor power Pauli group, which is 
generated through commutation by matrices of the form 
(Ta^VaK i,j £ {1, ... ,K}, a = x,y,z. 

The results discussed so far show that decoherence and 
decay can be completely suppressed in the limit of in- 
finitesimally short control time scale At — > 0. In order to 
convert this mathematical limit into a physically mean- 
ingful condition, we recall that complete information 
about the fluctuation-dissipation properties of a macro- 
scopic bath is encapsulated in the spectral density func- 
tion J{uj), measuring the density of modes at frequency 
lo multiplied by the square of the system-mode coupling 
strength. Quite generally, relaxation rates 7 = r~ e ] arise 
from an integration over the reservoir modes of the ef- 
fects due to thermal and vacuum fluctuations, weighted 
with J (to). This integration does not extend to arbitrar- 
ily large frequencies. For every physical spectral density 
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function, a finite ultraviolet cut-off frequency lo c always 
exists, such that J (to) — > for u> > lj c . The characteristic 
memory time r c ~ w" 1 , which is often set to be zero as a 
part of Markov approximation [jjj , determines the fastest 
time scale accessible to reservoir correlations. The con- 
dition At -C t c replaces, in a realistic scenario, the ideal 
limit At — > 0. Thus, the physical requirement for decou- 
pling is At < r c or lo c At < 1, i.e., the motion induced 
by the decoupling field needs to be faster than the fastest 
time scale characterizing the unwanted interactions. 

In the presence of a small but finite decoupling time 
scale At, the representation of the average Hamiltonian 

H in terms of the lowest-order contribution H ° is ap- 
proximate due to the higher-order terms. Consequently, 
decoupling is imperfect and relaxation dynamics still oc- 
curs even in the presence of control. For first-order de- 
coupling, the leading correction is due to and the 
key observation to estimate the decoupling accuracy as 
a function of At is that the overall coupling strength to 
the bath has been renormalized by a factor of the order 
uj c At. Since the spectral density depends quadratically 
upon the interaction strength, all the relevant relaxation 
effects tend to be suppressed by a factor of the order 
(u c At) 2 . For a generic fcth-order decoupling scheme, the 
controlled relaxation rate is able to be reduced as 



resources are still a stringent practical requirement, de- 
coupling techniques could be valuable compared to con- 
ventional error-correction networks in the field of NMR, 
ion-trap or solid state quantum computation. 

In summary, we showed how to manipulate the irre- 
versible component of open-system evolutions through 
the application of external controllable interactions. 
Maximal and selective decoupling were introduced within 
a common group-theoretical framework and their rele- 
vance to the issue of designing controlled effective open- 
system evolutions elucidated. In the spirit of weaken- 
ing the decoupling requirements as much as possible, the 
main question raised by the present analysis concerns the 
characterization of fault-tolerant decoupling schemes or, 
equivalently, the degree of decoupling attainable in the 
presence of power-limited and imperfect control opera- 
tions. Work in ongoing along these directions. 

This work was supported by ONR, AFOSR, and 
DARPA/ARO under the Quantum Information and 
Computation initiative and the NMR Quantum Comput- 
ing initiative. E. K. received support from the NSA. 
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(u c At) 



2 k 



At 

Tc. 



21; 



k > 1 



(13) 



Second-order decoupling can be realized via so-called 
symmetric cycles, whereby Ui(T c — t) — U\(t). Being 

H = for r odd, the error is improved to 0(A£/t c ) . 
Iterative pulse sequences for fcth-order decoupling can be 
designed for specific systems like quantum registers. 

In practice, the feasibility of the approach depends 
on both the relevant correlation times and the sophis- 
tication of the technology available to manipulate the 
specific physical system. In pulsed-NMR experiments 
p[ , where relaxation mechanisms due to "slow" nu- 
clear motions may involve correlation times longer than 
10~ 8 s, the main current limitation is represented by the 
pulse duration, rp « 1/is. In atomic physics, decou- 
pling methods could prove to be viable for damped har- 
monic oscillators schematizing the vibrational motion of 
trapped ions, since relevant cut-off frequencies may be 
estimated around 100 MHz and a variety of experimen- 
tal techniques exists for coherent optical manipulation 
As another potential area of applications, we men- 
tion semiconductor-based structures. Here, correlation 
times around Wp ei , se w 10 -13 s are comparable to the 
sub-ps time scale where control operations have been 
demonstrated Q and longer than the femtosecond scale 
of modern ultrafast laser technology. Rapid advance- 
ments in the capabilities of coherent control gives hope 
that, if not within the reach of present technology, im- 
plementations of decoupling schemes can be envisaged in 
a close future. In particular, since quantum computing 
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